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Abstract. Recurrence behavior of wave packets in coupled higher dimensional systems and periodically 
driven systems is analyzed, which takes place in the realm of higher coupling/modulation strength. We 
analyze the wave packet dynamics close to nonlinear resonances developed in the systems and provide the 
analytical understanding of recurrence times. We apply these analytical results to investigate the recurrence 
times of matter waves in optical lattice in the presence of external periodic forcing. The obtained analytical 
results can experimentally be observed using currently available experimental setups. 
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1 Introduction 

In one dimensional quantum systems with discrete en- 
ergy spectrum, quantum recurrence theorem establishes 
the existence of recurrent dynamics [T]. Quantum char- 
acteristics of dynamical systems, which exhibit chaos in 
their classical domain, have posed interesting questions 
for researchers. These systems in general display inher- 
ent quantum interference phenomena which modify the 
quantum dynamics and produces a contrast in their evo- 
lution. Quantum interference phenomena, constructive or 
destructive in nature, contribute to the quantum recur- 
rences in dynamical systems. Dynamical recurrences in 
such systems originate from the simultaneous excitation 
of discrete quasi-energy states [5]. These recurrences are 
suitable to use as a probe to study quantum chaos [3)111 
[5]. In this contribution we extend the discussion on the 
quantum recurrence phenomena and treat them in higher 
coupling/strong modulation regimes and study the quan- 
tum dynamics of the system for the nonlinear resonances. 
Our theoretical work is based on Floquet theory , which 
presents an elegant formalism for the analysis of peri- 
odically driven and coupled higher dimensional systems 
[?]■ Floquet state formalism has been applied to a num- 
ber of time-dependent problems: from coherent states of 
driven Rydberg atom |8j, chaotic quantum ratchets 
electron transmission in semiconductor hetero-structures 
[To], selectively suppressing of tunneling in quantum-dot 
array [TT] to frequency-comb laser fields jl2 . The Floquet 
states have also been used in Bose-Einstein condensate 
systems and applied to probe superfluid-insulator tran- 
sition [13], towards coherent control [TJ] and dynamical 
tunneling [IS] . 
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In this paper: (i) We extend Floquet analysis for non- 
linear resonances and explain the existence of robust dy- 
namical recurrences in dynamical system, in contrast to 
delicate dynamical recurrences already reported [SlI^fTBl 
I17j ; (ii) Super revival times in dynamical systems are dis- 
cussed first time for both delicate dynamical recurrences 
and robust dynamical recurrences; (iii) We show that, the 
non-linearity of the uncoupled systems, and the initial con- 
ditions on the excitation contribute to the classical period, 
quantum revival time and super revival time occurring 
in the coupled higher-dimensional systems or periodically 
driven systems; (iv) We apply these results to the dynam- 
ics of matter waves in driven optical lattices, a topic of 
current research [TSIIT^I^ . Classical dynamics of the sys- 
tem display an intricate dominant regular and dominant 
stochastic dynamics, one after the other, as a function of 
increasing modulation amplitude [^122] . 

The paper is organized as follows: The effective Hamil- 
tonians for nonlinear resonances in higher dimensional sys- 
tems (Sec. 12.11) and periodically driven time dependent 
systems (j2.2p are derived. In Sec. [31 we explain the time 
scales and their dependence on system parameters. We ap- 
ply obtained results to the matter wave dynamics in driven 
optical lattices and explain the time scales as a function of 
system parameters. We compare analytical and numerical 
results in Sec. |H 



2 Floquet Theory of Nonlinear Resonances 

Floquet theory is an elegant formalism to quantize coupled 
higher dimensional systems (HDS) and periodically driven 
system (PDS) [7]. In HDS, with two degrees of freedoms 
coupled, such as Billiards [53] and in PDS, we observe var- 
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ious dynamical modes. In these Hamiltonian dynamical 
systems the stable nonlinear resonances are immersed in 
stochastic sea and the system may display global stochas- 
ticity beyond a critical value of coupling or modulation 
strength pm^SJUS]. Here, we apply Floquet analysis to 
find quasi energy eigen states and quasi energy eigen val- 
ues in HDS and PDS and discuss their parametric depen- 
dencies on the uncoupled system and coupling constant or 
modulation strength. 



2.1 Higher Dimensional Systems 

We write general Hamiltonian for a system, with its de- 
grees of freedom coupled, that is. 



and 



4AF 



(5) 



The parameter C, defines the nonlinear dependence of the 
energy of the system on the quantum number and S is 
re-scaled Plank's constant. 

The TT-periodic solutions of equation ^ correspond 
to even functions of the Mathieu equation whose corre- 
sponding eigenvalues are real [33 . These solutions are 



defined by Floquet states. 



e'l^'^Pf.ie), where. 



tt), and fi is the characteristic exponent. 



H = Ho{I) + \H,{I, 



(1) 



where, Hq is the Hamiltonian of the system in the ab- 
sence of coupling, expressed in the action coordinates / — 
(/i,/2). Moreover, He is the coupling Hamiltonian which 
describes coupling and is periodic in angle, — {01,62) 
and the parameter A is coupling strength. The Hamilto- 
nian may express Billiard systems in general systems 
and same time is a good candidate to describe multi- 
atomic molecule |^ and Lorentz gas Jl^ . We express 
the coupling Hamiltonian as, 



In order to write the 7r-periodic solutions in (/5-coordinate, 
we require /i to be defined as /i = ^{j) = 2j/N, where, 
j=0,l,2,....,iV-l. 

The allowed values of can exist as a character- 
istic exponent of solution to the Mathieu equation for 
discrete i> (which takes integer values) only for certain 
value a^{pij),q), when q is fixed. Hence, with the help of 
Eq. (|4]), we obtain the values of unknown £. Therefore, we 
may express the eigen energy of the system as [55] 



8 



mod Huj. 



(6) 



Hnil)e 



in. 6 



where, n = (ni, 712). Whenever, the frequencies J7 ■ 



Here, a = M/N defines the winding number. 
/2') In order to check this result we study the case of zero 
coupling strength, that is A = 0. In this case, the value for 
Mathieu characteristic parameter becomes a,^{q = 0) = 
dHo/dl rpj^jg reduces the quasi-energy f^.i/, in the absence of 



obey the relation, n.H = 712^^2 = 0, resonances oc- coupling term, that is, A = 0, provided z/ as 2(/-/o)/iVS. 

cur in the system. 

We consider that the coupling exists between actions 
Ii and l2- Within the region of resonance, we find slow 
variations in action, hence, following the method of secu- 
lar perturbation theory, we average over faster frequency, 



2.2 Time Dependent Systems 



and get the averaged Hamiltonian for the N resonance, 
asH = Ho{I) + XVcos{Nip_). Here, ip^ = ei-{M/N)02, 
/ = /i is the action corresponding to the angle 9i, V is the 
Fourier amplitude, and M and TV are relatively prime in- 
tegers |24| . Averaging over rapidly changing ip^ makes the 
corresponding action I2 as the constant of motion. More- 
over, Ho{I) expresses uncoupled averaged Hamiltonian. 



In the periodically driven potentials energy is no more a 
constant of motion. For the reason we solve the time de- 
pendent Schrodinger equation by using secular perturba- 
tion approximation as suggested by Max Born j36l . There- 
fore, the solution is obtained by averaging over rapidly 
changing variables. This leads us to find out a partial so- 
lution of the periodically driven systems. As a result, we 
find quasi energy eigen-states and the quasi eigen energies 



The energy of initially narrowly peaked excitation change^f the dynamical system for non-linear resonances. 



slowly when we produce it in the vicinity of TV*'' non-linear 
resonance. Therefore, we expand the unperturbed energy, 
Hq (/) , by means of Taylor expansion around mean action 
I — Iq and keep only the terms up to second order 

We [5U1I5T] use operator definition of action, that is, 

I = j-gg to get the near resonant energy spectrum of the 
system by mapping the corresponding Schrodinger equa- 
tion on the Mathieu equation [32 , where, the new variable 
= Nlp/2, hence, we get 



In order to study the quantum nonlinear resonance of 
the TDS, we consider the Hamiltonian, 



H = Ha + XV{x) sm{t). 



(7) 



The solution of the Schrodinger equation corresponding 
to the Hamiltonian ([7]) can be written in the form [ro[|29[ 



01 



a-2qcos2e](l>{0) =0, (3) 

where, the Mathieu characteristic parameter is 

8 - ,.,2 



n) 



t t 
7V% 



}• (8) 



2C 



(4) 



Here, En is the mean energy, C„(t) is time dependent 
probability amplitude, n is mean quantum number and 
\n) are eigen states of undriven system. On substituting 
equation ([5]) in the time dependent Schrodinger equation. 
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we find that the probabihty amplitude C„(i), changes with 
time following the equation, ikCn{t) — [En ~ En — (n — 
' ]Cn[t) + ^[Cn+N - Cn-N). Where, V = Vn 



-N 



Vn+N are off-diagonal matrix elements and are approxi- 
mately constant near the potential minima for tight-binding 
approximations . 

We take the initial excitation, narrowly peaked around 
the mean value, n. For the reason, we take slow variations 
in the energy, Em around the fi in a nonlinear resonance, 
and expand it up to second order in Taylor expansion. 
Thus, the Schrodinger equation for the probability ampli- 
tudes, Cn{t), is 



n){uj~j^)Cn{t) + ^T:{n- 



n)\Cn{t) 



Cn 



In order to obtain equation ([9]), we average over the fast 
oscillating terms and_^keep only the resonant ones. Here, 

and the non-linearity, in the 



the frequency, uj = ^ 



dE„ 

itdn \n=n-, 



time independent system is given as C = Ji^"^ \n=n- 



We introduce the Fourier representation for C„(i) as, 
Cn{t) = T7^ g(0,t)e-'("-'"')«/^d0, (10) 



2N7r J, 

which helps us to express equation ^ as the Schrodinger 
equation for g{6, t), such that ikg{6, t) — H{9)g{9, t). Here, 

the Hamihonian H{0) is given as, H{e) -^^^S^ - 
iNk{u: — jf)-§g — Ay sin 6*. In order to obtain this equa- 
tion, we consider the function g{9,t) as 2Ntt periodic, in 
9 coordinate. 

Due to the time-independent behavior of the Hamil- 
tonian, we write the time evolution oi g{9,t), as g{9,t) = 

~ —let 

g{9)e . Therefore, Schrodinger equation for g{9,t) re- 
duces to the standard Mathieu equation (jS]). Here, we 
express g — x(z) exp (— i2(A^aj — l)z/iV^C^), where, 9 — 
2z + 7r/2. The same as in equation ([3]), the Mathieu char- 
acteristic parameter, a, is defined as 



8 ANuj-lf 



and 



(11) 



(12) 



Hence, the quasi eigen energy of the system, obtained from 
equation ([TT|). is 



f — 



mod JeL 



(13) 



where, the index ly takes the definition, v = ^''"^y , where, 
n is the quantum number for time independent system. 

3 Quantum Revival Times Based on Quasi 
Energy Spectrum 

In this section we discuss two cases correspond to non- 
linear resonances. We consider weakly coupled g < 1 or 



strongly coupled g 3> 1 potentials. In g < 1 situation, for 
large v |16j and in q 3> 1 situation for small i^, near the 
center of resonance matrix elements are constant. In the 
following discussion, we analyze the wave packet dynamics 
in these regimes and show their parametric dependencies. 

The time scales, T^^^ at which recurrences of an ini- 
tially well-localized wave packet occur depend on the quasi- 
energy of the respective system. The recurrence times are 
obtained as, T*^^^ = 2n/Q^^\ hence, the values of j as 
j = 1,2,3..., correspond to, respectively, classical, quan- 
tum, super, and higher order revival times. With the help 
of Eqs. © and (fT5)) . we obtain the frequencies J7*^^^ as 



^(1) ^ j_ (d£^ ^ d£^^u 
dv 



(9) = 



^ 1 f (^^^i-iM ^ ^ d^^tJ.,v ^ ^ d'^£i-LM ^ d^£^,^ 
After a few mathematical steps Eqs. (ITil) reduce to 



(14) 



I (d^£., , ^2d£.. 



2\k^ 5i/2 



ak- 



du 



(15) 



which lead to classical period, T*-^' — , quantum re 
vival time, T^^) = and super revival time, T'^'^' = 

rjn{spr) ^ when calculated at the mean values. 



3.0.1 Delicate Dynamical Recurrences 

The condition, q < 1, may be satisfied in the presence of 
weak perturbation due to external periodic force (and/or), 
for large nonlinearity and/or for large effective Plank's 
constant 16,5,2,17,. The Mathieu characteristic parame- 
ters, fljy and bv are given [33], as 



2(1/2 _ 1) 



for 1/^5 



(16) 



The above expressions arc not limited to integral value 
of V and are very good approximations when v is of the 
form, m -|- ^. In case of integral value of i.e. u — m, the 
series holds only up to the terms not involving — in 
the denominator. 

The energy spectrum for weakly modulated periodic 
potentials [S^ can be defined using equations ([T3| and 
p^ . The relations obtained for classical period, quantum 
revival time and super revival time, in the presence of 
small perturbation for primary resonance N — 1, index j 
takes the value, j — 0, and time scale are simplified as 



(cl) 



(1 



(17) 
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j^irev) ^^^_^^^rev)^r^{rev)^ (18) 
(spr) _ TTO;^ (1 - )4 

~ 2XVCA^ /.I ' ^^^^ 

here, T^"'^ = 77 classical period and tI^"'"'''' = 27r/(ifeC) 
is quantum revival time in the absence of external mod- 
ulation. Furthermore, A = (1 — —). The modification 
factors are given as 

2^ c^2 ^ 



For primary resonance iV = 1, index j takes the value, 
j = 0, and time scale are simplified as 



id) 



(c/) 



8^ 



(4 + /ii)v/C (4 + m)\ 



-7(rei!) 



(re-t;) 



1 - 



(8/^l)2g 



3(mi+4) , 9(4 + ^i)2 



and 



and 



, (25) 



(26) 



(27) 



2' W2 ^ (1-^2)3' 

here, /ii = is re-scaled non-linearity. 
3.0.2 Robust Dynamical Recurrences 

The condition, g 3> 1, may be satisfied in the presence 
of large amplitude. A, of external modulation by periodic 
force in a dynamical system. In addition, we may get the 
regime by considering a system with very small linearity 
i.e. C ~ and/or by taking effective Plank's constant be 
very small. Quasi energies for nonlinear resonance are de- 
fined in terms of Mathieu characteristic parameters, which 
are 



-2q + 2sVg - 



1 s-^ + 3s 



23 



(20) 



where, s = 2v + 1. 

The energy spectrum of nonlinear resonances for strongly 
modulated periodic potentials can be defined using equa- 
tions (IT^ and (PO)) . Here, in the deep potential limit, the 
band width is. 



24..+5^/|^f + 1 cxp(-4^) 



+1 



(21) 



Atoms in modulated optical lattice: In order to analyze the 
time scales in equations (IT7|) - (ITO1) and equations (j25ll27p . 
we consider ultracold Rubidium atoms, Rb^"^ , in standing 
wave field with phase modulation due to acousto-optic 
modulator. The system is of immense theoretical and ex- 
perimental interest since last two decades in atom optics 
with cold atoms [22 and Bose-Einstein condensate |38) . 
The atoms in the presence of phase modulation experience 
an external force, thus exhibit dispersion both in classical 
and quantum domain. In quantum dynamics atoms expe- 
rience an additional control due to effective Plank's con- 
stant. This as reported earlier limits the classical evolution 
observed in classical prototype and results in dynamical 
localization [39]. Here, we report that in long time evo- 
lution, material wave packet display quantum recurrence 
phenomena. The dynamics of cold atoms in phase modu- 
lated optical lattice is governed by the Hamiltonian, which 
is 

H = ^ + Y cos[2fcL{a; - ALsm{ujjnt)}], 

where, fc/, is wave number and Vq define the potential 
depth of an optical lattice. Furthermore, AL and uJm are 
amplitude and frequency, whereas, M is the mass of an 
atom moving in the optical lattice. 

The unitary transformation^ 
4>{x,t) = Vi(x, t) exp[i{ ""^^^^^^°^("-*) X + m}l where. 



Keeping lower order terms in s, in equation (j20[) . we 
get quasi energy spectrum for nonlinear resonances as 
■^)'kuJh, which resembles to the harmonic oscillator energy 
spectrum for ujh = 2-\/Vo- 

Similarly the relations for classical period, quantum 
revival time and super revival time for strongly driven 
case are obtained as, 



the lattice, modified the Hamiltonian as 



/?(i) ~ "m^^ ^" |- sin^zi^„r; ^ ^ frame co-moving with 



Fx sinujmt, 



(28) 



T. 



(cl) 



An 



(22) 



{rev) 



327r 



^'C{l + ^^|^+8a^Vg-(2f + 1) ak} 



(23) 

j^jspr) _ 327r ^, [11 3"fe(2^' + l) ij^ (24) 



where, F = MALoj'^ is amplitude of inertial force emerg- 
ing in the oscillating frame. In laboratory accessible quan- 
tities, scaled optical lattice depth ^ is 5 to 20. For exam- 
ple Rb^'^ atoms trapped by optical lattice of wavelength 
842nm has Er — 1.34xl0~^^ey, thus recoil energies are 
of the order of 10~^"eV^ 40 . To examine the dynamics of 
cold atoms in driven optical lattices, Hamiltonian ((28)) is 
expressed in dimensionless quantities. Using scaling trans- 
formations z = k^x, T — uJmt, ^{z, t) — ^{x, t) and mul- 
tiplying the Schrodinger wave equation by where, 

^ The unitary, time periodic transformations preserves the 
quasi-energy spectrum. 
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LUr = ym ' single photon recoil frequency, we get dimen- 
sionless Hamiltonian, 



2 9z2 



— cos Zz + \z sm T, 



here, A = kL^L, is modulation amplitude, i.e. the effective 
amplitude of inertial force, Vq = and r is scaled time 
in the units of modulation frequency utm- In this case the 
rescaled Plank's constant k — 

As time scales of driven optical lattice are expressed 
in terms of undriven system parameters, for the reason, 
we investigate, classical period, quantum revival time and 
super revival time scales of undriven optical lattice [ITll^ 

sg. 



In the case of undriven shallow optical potential Tq'^^^ = 

(1 + 5(Af)^' ^o"^"^ - 2.(1 - 4{St?) and T^-^^ = 
^^},\ , where, n is band index of undriven lattice, while, 
go = is rescaled potential depth in units of single re- 



2M 



coil energy 

cal lattice, the classical time period is Tq"' = 



On the other hand for deep opti- 



3(s" + l) 



syfj ' 2** go -^^ where, s — 2n + 1. The quantum re- 
vival time xlf^^'^ = 47r(l — -^^) and super revival time is 



327r^ 



rp(spr) 

In shallow lattice potential limit, i.e. qo < 1, neglecting 
the higher order terms in qo the classical frequency, uj = 

2 

'2,n{l ~ „, and non-linearity, C = 2 



% 3ra^ + l 
2 (n^-l)i' • 



In deep optical potential limiting case, i.e. go 3> 1, the 
classical frequency is a; = 4(.ygo ~ ^"g^"*" ) and non-linearity 



3(2rt+l) I 



is C = I - 1 - 

For driven optical potential, keeping in view the pe- 
riodicity of Floquet solutions [121111], we take only even 
values of the index [33] . Hence, re-scaled Mathieu charac- 
teristic exponent, v = 2{l + /?), where, (3 = ^"^^ . 

Comparing the coefficients of eigen energy of the un- 
driven system and equation of motion of probability am- 
plitude in the absence of modulation, we get / = , 
which is new band index for nonlinear resonance and at 
the center of resonance I — 0. 

For delicate dynamical recurrences, time scales for pri- 
mary resonance iV = 1, are 



rpicl) rTi{cl) 



q 



1 



2 {4(? + /3)2-l}2 



(rev) rp{rev) 



T. 



and t['^''^ = 



1 - 



12G + /3)2 



1 



2 {4(/ + /?)2 - 1}3 

^{4(; + /3)2„i}4 



iN'^Qkq^il + + PY + 1} ■ 



(29) 

(30) 
(31) 



Where, Tr 



icl) 







[rev) _ 4^ 



is the classical time period and 



T' 
tern. 



^ is quantum revival time for unmodulated sys- 



On the other hand, in deep optical lattice limiting case, 
time scales for the atomic wave packet for primary reso- 
nance with A'^ = 1, are given as 



p(cO _ 



2tt 



-^{rev) 



N^kCiVq- 



Stt 



4{1 + I3) + 1 



} 



N^kC 



and 



T. 



(spr) _ 327r y/q 
mkC ' 



(32) 



(33) 



(34) 



In case of deep lattice, when external modulation fre- 
quency is close to the harmonic frequency, matrix ele- 
ments, V can be approximated by those of harmonic os- 
cillator and Mathieu parameter q can be approximated as 
4Vn+iA rr^ Under this approximation time scales 



q 



are 



rp{clj _ 



167:q^/N^VC 



16(n + l)iVX- {4(Z + /3) + Ijq^k^ 



(35) 



(rev) 



fPkC 



and T 



(spr) 
A 



3{4(/ + /3) + l}go^feVC 
647r(n-h 1)3%/A 



(36) 



Behavior of classical periods, quantum revivals and su- 
per revivals of matter waves in modulated optical crys- 
tal in nonlinear resonances versus modulation is shown in 
Fig-[T] In each plot of this figure, left vertical axis shows 
the time scales when shallow optical lattice is weakly or 
strongly modulated and right axis shows the time scales 
when deep optical lattice is weakly or strongly modulated. 
The upper row of Fig-[1] represents the time scales for small 
q values i.e. delicate dynamical recurrences as a function of 
modulation A, while, lower row represents the time scales 
when g ^ 1, as a function of modulation A, i.e. robust dy- 
namical recurrences. Here, left column shows the results 
related to the classical periods. Quantum revival times are 
plotted in middle column, while, right column shows super 
revival times. 

We note that when optical lattice is perturbed by weak 
periodic force, the classical period increases with modula- 
tion, as given in equation ([^ . Classical period for weakly 
driven shallow lattice potential changes slowly as com- 
pared to weakly driven deep lattice potential as shown in 
Fig-lIJa) . When lattice is strongly modulated by an exter- 
nal periodic force, the classical period decreases as modu- 
lation increases. Classical period for strongly driven opti- 
cal lattice is given by equation ([5^ . The behavior of clas- 
sical period for strongly driven lattice versus modulation 
is qualitatively of the same order for both strongly driven 
shallow lattice and strongly driven deep lattice as shown 
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in Fig-[IJb). The behavior of classical period in strongly 
driven lattice case is understandable as strong modulation 
influence more energy bands of undriven lattice to follow 
the external frequency and near the center of nonlinear 
resonance the energy spectrum is almost linear, as can be 
inferred from equations and (PH]) with assumptions 
q ^ 1 and / is small. 

Quantum revival time in nonlinear resonances versus 
modulation is shown in middle column of Fig-[TJ For deli- 
cate dynamical recurrences Fig-[TJc), the quantum revival 
time decreases as modulation increases. The behavior of 
quantum revival time is given by equation (1301) and equa- 
tion p3p respectively for delicate dynamical recurrence 
and robust dynamical recurrence. For weakly driven shal- 
low optical lattice or weakly driven deep lattice, qual- 
itative and quantitative behavior of revival time is al- 
most similar. The qualitative behavior of revival time for 
strongly driven shallow lattice is different from the strongly 
driven deep lattice Fig-[TId), as in the later case, change 
in revival time is almost one order of magnitude larger 
than the former case, for equal changes in modulation. 
Classical period and quantum revival time for delicate dy- 
namical recurrences show good numerical and analytical 
resemblance for the system with our previous work [5 ,1?1 
[32] ■ Here, the difference in the revival time behavior for 
strongly driven shallow lattice as compared to strongly 
driven deep lattice is due to the difference in energy spec- 
trum of undriven system. In deep lattice, due to small non- 
linearity more energy bands are influenced by the external 
drive and resonance spectrum is similar to that of har- 
monic oscillator near the center of nonlinear resonance. As 
modulation is increased more and more energy bands are 
influenced by external drive. Fig-[5] shows spatio-temporal 
evolution of an initially well localized wave packet in a 
lattice minima. Upper panel is for the spatio-temporal dy- 
namics of atomic wave packet in the absence of periodic 
modulation, while the lower panel presents the case when 
external modulation, X ~ 3, q — 85.14. Time evolution 
of wave packet in optical lattice shows that wave packet 
splits into small wavelets and spreads over the neighboring 
lattice sites. Later, these wavelets constructively interfere 
and wave packet revival takes place. It is clear that with 
the introduction of modulation, revival time changes as 
shown in Fig-[T] which is due to a change in interference 
pattern and confirms analytical results as discussed above. 

The super revival time behavior versus modulation, 
shown in right column of Fig-[l] for delicate dynamical re- 
currences and robust dynamical recurrences are given re- 
spectively in equations (|5T|) and (IMl) . As lattice is slightly 
perturbed by periodic external force, the super revival 
time, T^*'"^-' decreases with modulation for delicate dy- 
namical recurrences as shown in Fig-[IJe). On the other 
hand, the super revival for robust dynamical recurrences 
increases with modulation as shown in the Fig-[IJf). Here, 
qualitative behavior of super revival time is same but 
quantitatively super revival time increases almost two times 
faster. 



4 Discussion 

We presented general discussion on the occurrence of ro- 
bust dynamical recurrences in higher dimensional and time 
periodic systems in the presence of strong coupling/modulation. 
The dependence on energy spectrum of time scales is ex- 
plained. We applied our results to the atomic dynamics in 
optical lattice driven by the periodic forcing. 

We have two cases for each condition i.e. for q < 1 
condition which is satisfied either shallow or deep poten- 
tial is slightly perturbed by small external periodic force. 
Analytical expressions (|17m9p for wave packet time scales 
are valid. Parameters, V, to and C are scaled matrix ele- 
ments, frequency and non-linearity of respective undriven 
shallow or deep potential. Similarly, condition g 3> 1 is 
satisfied when shallow or deep potential is strongly mod- 
ulated. Expressions (I25II27P represent time scales of the 
driven system in this case. Moreover, for deep lattice case 
(undriven or driven), expansion (j20|) is quite good for 
Mathieu characteristic exponent i' satisfying the condition 
ly ^ and for i/ ^ ^/q, equation (fTC)) , satisfy the numer- 
ically obtained energy spectrum. The intermediate range 
(i^ ~ ^/q), where, energy spectrum changes its character 
from lower to high value is estimated as ~ 2 || -^/l || 
[45] . here, w y w denotes the closest integer to y. 

We numerically explore the dynamics of quantum par- 
ticle in a modulated optical lattice for a nonlinear res- 
onance. Numerical results are obtained by placing wave 
packet in a primary resonance with iV = 1 . This is realized 
in experimental setup of Mark Raizen at Austin, Texas. 
The authors [3^ worked with sodium atoms to observe 
the quantum mechanical suppression of classical diffusive 
motion and employed the {3Si,F = 2) — > (3P3,F = 3) 

transition at 589nm , with a;o/27r = 5.09 x lO^^Hz. The 
detuning was Sl/^tt = 5.4 x 10^ Hz. The recoil frequency 
of sodium atoms was a;r/27r = 25kHz for selected laser 
frequency and modulation frequency was chosen a;m/27r = 
l.iMHz, whereas, the other parameters were k = 0.038 
and g = 55 (or Vo — 0.16). In the other experiment 
|46| . Bloch oscillations of ultracold atoms were observed 

with the 65*1 6P3 transition in cesium atoms setting 

22 

Al = 852nm and a;o/27r = 3.52 x 10"i7z. The detuning 
was Sl/^tt = 3 X lO^^Hz, with q up to 1.5. The recoil 
frequency was ujr/2Tr — 2.07kHz, so that a driving fre- 
quency uim/2iT = lO^Hz, three orders of magnitude lower 
than in the former experiment, gives S s:! 4, taking the 
dynamics to the deep quantum regime. We numerically 
investigate the validity of our results with known theoret- 
ical [21] and experimental [39 results with dimensionless 
rescaled Plank's constant k = 0.16, Vq = 0.36, which is a 
case of strong modulation to a deep optical lattice of po- 
tential depth Vq — 28.13Er. The wave packet is initially 
well localized in such a way that localization length is less 
than or order of lattice spacing. 

To have an idea about the classical dynamics of the 
system, we plot the Poincare surface of section for mod- 
ulation strengths A = 0, 1.5 and 3.14 as shown in Fig-[3l 
From phase space plot, we see the appearance of 1:1 res- 
onance for A > 0. This resonance emerges when the time 
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period of external force matches with period of unper- 
turbed system. One effect of an external modulation is 
the development of stochastic region near the separatrix. 
As modulation is increased, while the frequency is fixed, 
the size of stochastic region increases at the cost of regular 
region. 

In order to observe the dynamics of a quantum parti- 
cle inside a resonance, we evolve a well localized Gaussian 
wave packet in the driven optical lattice. Square of auto- 
correlation function (| C |^) for the minimum uncertainty 
wave packet is plotted as a function of time in Fig-|31 for 
A = 1.5,(7 = 27.78, whereas, Ap = Az = 0.5, = 0.5 
and Vq = IQEr- In Fig-|3] upper panel is an enlarged view 
which displays classical periods, while middle panel of the 
figure shows quantum revivals. Lower panel shows the ex- 
istence of super revivals. The classical period, quantum 
revival time and super revival time are indicated by ar- 
rows and showing there characteristics. Numerical results 
are in good agreement with analytical expressions. 
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Fig. 1. Left panel: Classical time period versus A for weak external modulation (a) and for strong external modulation (b). In 
this figure, for deep lattice Vb = 16-Br-, for shallow lattice Vo ~ 2Er and S = 0.5. Middle panel; Quantum revival time versus 
A for weak external modulation (c) and for strong external modulation (d). Right panel: Super revival time versus A when 
external modulation is weak (e) and for strong external modulation (f). 
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Fig. 2. Spatiotemporal dynamics of atomic wave packet in the absence of modulation (a), and in the presence of modulation 
(b). The amplitude of modulation is A = 3 and q = 85.14. Other parameters are Vb = 0.36, Ap = 0.1 and H = 0.16. 




Fig. 3. Poincare phase space for driven optical lattice for modulation strengths A = 0, 1.5, 3 and potential depth Vb = 16Er 




Fig. 4. Square of auto-correlation function of a Gaussian wave packet is plotted versus time. Here, the parametric values are 
A — 1.5, — 0.5, Vq = 16Er, Az = Ap = 0.5. The wave packet is initially well localized at the central lattice well around the 
second band of undriven lattice. Classical period, revival time and super revival time are seen and indicated by arrows. 



